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LIGHT REFLECTION IS NONLINEAR OPTIMIZATION 



JIAKUN LIU 



Abstract. In this paper, we show that the near field reflector problem is a nonlinear 
optimization problem. From the corresponding functional and constraint function, we 
derive the Monge- Ampere type equation for such a problem. 



1. Introduction 

Qh ■ Optimal transportation, due to its various applications, has been extensively studied 

in recent years. The modern theory of optimal transportation is mainly built upon Kan- 

<£ torovich's dual functional, which is a linear functional subject to a linear constraint. With 



his dual functional, Kantorovich introduced linear programming, which is a class of linear 
optimization problems. An important new application is the reflector design problem. In 
], Xu-Jia Wang showed that the far field case of the reflector design problem is an optimal 



ctf 



transportation problem, and so is a linear optimization problem. The purpose of this paper 

O ■ 

l/-j ■ is to show that the general case of the reflector problem is a nonlinear optimization prob- 

cn 

en 



lem. More examples of nonlinear optimization problems and also questions of the existence 
and regularity of potential functions and optimal mappings will be investigated in [8] and 
subsequent papers. 

Suppose that a point source of light is centered at the origin O and for each X S C § n 
we issue a ray from O passing through X, which after reflection by a surface T will illuminate 
^ ! a point Y on the target surface ST in R n+1 . Let / e L l (£l),g E L 1 ^*) be the input and 

gain densities, and dfi, dv denote the surface area elements of £1, Q*, respectively. The near 
field reflector problem can be formulated as follows: given (£l,f) and (f2*,g) satisfying the 
energy conservation condition 



(1.1) fdfji= gdv, 

Jn Jn* 

find a reflector T such that the light emitting from 0, with density / is reflected off T to the 

target £1* and the density of reflected light is equal to g. 
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In our reflector problem, we assume that both fi and Q* are compact and each has 
boundary of measure zero. Represent the reflector T as a radial graph of function p, 

(1.2) r = {Xp(X) : X G 0}. 

Let V{p, v) be the set of measures on Q x Q* with p, v as their marginals. Let 7 E V(p, v). 
Denote by C+(f2) the set of positive continuous functions on O. Define a functional 

(1.3) I(u,v)= [ F(X,Y,u,v)d 1 , 
for (u,v) £ C+(Cl) x C+(tt*), where 

(1.4) F(X,Y,u,v) = f(X)logu + g(Y) (logv + log(l {X /J\ =)) 

and ( ,) is the inner product in M n+1 . The main result is the following: 

Theorem 1.1. Suppose that f,g are two bounded positive functions on f2,f2* ; respectively, 
such that (jl.ip is satisfied. Suppose that f2* is contained in the cone Cy = {tX : t>0,l£ 
V} for a domain V C S n and 

(1.5) HnF = 0, 

where Q and V denote the closures o/O and V , respectively. Then there is a dual maximizing 
pair (p, n) G K, which satisfies 

I(p,rj)= sup I{u,v), 

(u,v)£K 

where I(u,v) is given in (|1.3p - (|1.4p . and the constraint set K, is given by 

lC = {(u,v) G C + (n) xC + (C') : <p{X,Y,u,v) < 0} , 
with the constraint function 

(1.6) </>(X,Y, U ,v) =logu + logv + log(l ( ^ : ) . 

Moreover, p is a solution of the reflector problem with given densities (Q, f) and (Q*,g). 

In Theorem ll.il the functions p, n are also called potential functions, and a solution of the 
reflector problem needs to be understood as a weak solution. The notion of weak solutions 
was introduced in (6j[7], see §2.2 below. It follows from Remark 12.11 that for each choice of 
the parameter cq > 0, there is a weak solution p satisfying info p > cq. 

Moreover, we show that the function p solves a Monge-Ampere type equation. Assume 
that Q* is given implicitly by 

(1.7) n* = {Z € R n+1 : 4){Z) = 0}. 
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Suppose that O is a subset of upper unit sphere S™ = S n n {x n+ i > 0}. Let X = (x, x n+ \) 
be a parameterization of 0, where x n+ \ = yl — |x| 2 =: u>(x), and x = (xi, ■ ■ ■ ,x n ). For 
simplification, we define some auxiliary functions 

(1.8) a = \Dp\ 2 -(p + Dp-x) 2 , 

(1.9) b = \Dp\ 2 + p 2 -(Dp-x) 2 , 

(i.io) t = p x ^~y^\ (3 



pXn+l ~ (Y-Xp)-Vil>' 

and denote the matrix 

(1.11) M = {M l] ], M tJ =8 l3 + 



% ' 1 19' 

1 — \x\ z 
By computing in this local orthonormal frame, we obtain our equation as follows 

Theorem 1.2. The function p is a solution of 



(1.12) 



det 



D 2 l D ® D P _<}_A M 

p 2tp 



a n+1 



t n bf3 



f 



2 n p 2n + 1 u 2 g\V^\' 



The equation (|1.12p was previously obtained by Karakhanyan and Wang studying the 
near field reflector problem [6]. One of the main differences in our derivation of (|1 . 12[) is 
that instead of applying the reflection law as in [6], we have differentiated the constraint 
function (jl.6p directly in general cases, see (|3.4p below. We remark that our method is more 
general and can be applied to the study of other reflector and refractor problems, |8]. 

This paper is arranged as follows. In Section 2, we first introduce a class of nonlinear 
optimization with potential functions, and then prove Theorem II .11 In Section 3, we derive 
the equation for potentials arising in general nonlinear optimization problems, and apply 
this formula to prove Theorem 11.21 In Remark 13.11 we point out that the far field reflector 
is a limit case of the near field one and related to a linear optimization problem. 

2. Formulation to optimization 
2.1. Nonlinear optimization. In general, we consider a functional 



(2.1) I(u,v) = F(x,y,u,v)dj, 

JUxV 

for (u,v) & C(U) x C(V), where U,V are two compact domains in W 1 or a manifold Ai n , 
F is a function on U x V x IR 2 , and the measure drj has marginals dx,dy, which are the 
volume elements of U,V, respectively. 

We want to maximize the functional / among all pairs (u, v) in a constraint set 

(2.2) K = {(u,v)e C{U) x C(V) : <p(x, y, u,v)<0 in U x V} , 
where d> is the constraint function defined on U x V x M 2 . 
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Note that when 

(2.3) F(x, y, u, v) = — f(x)u(x) + -r^rg(y)v(y), d-y = dxdy, 

for some / > 0, G ^(U), g > 0, G i 1 (V A ) satisfying f u f = f v g, and 

(2.4) 0(x, y, u, u) = u(x) + v(y) - c(x, y), 

we have a linear optimization problem related to an optimal transportation with the cost 
function c, and mass densities f,g supported on U, V, respectively. See [T| |2| [TT| [T2] . 

Definition 2.1. A pair (u,v) G tC is called dual pair with respect to 4>, if 

(2.5) u(x) = sup{t : 4>{x,y,t,v{y)) < 0, Vy G F}, 

■u(y) = sup{s : 4>(x , y , u(x) , s) < 0, Vx G [/}. 

If furthermore I(u,v) = sup^I, (u,v) is called dual maximizing pair of I. In such a case, 
u, v are also called potential functions in the nonlinear optimization (|2.ip — (| 2 . 2 j) . 



Write F = F(x,y,t,s) and <fi = <p(x,y,t,s), where x,y,t,s are independent variables. 
Use the subscripts to denote the partial derivatives, i.e. F t = dF/dt, cj) s = d<f>/ds, etc. We 
always assume that F is C 1 smooth in t, s and integrable in x, y; 4> is C 2 smooth in all 
variables. Moreover, we assume the following conditions on F and <f>: 

(i) F(x,y,t,s) is monotone increasing in t, s, namely 

(2.6) F t >0, F s >0, y(x,y,t,s) G Ux V xMxi 

(ii) 4>(x, y, t, s) is strictly increasing in t, s, namely for a constant 5$ > 

(2.7) (f>t>S Q , (f> s >6 , V(x,y,t,s) G U x V xlxl. 
(iii) for any pair (u, v) G fC, the balance condition holds: 



(2.8) / 1-F t (x,y,u(x),v(y)) + F s ^(x,y,u(x),v(y))\d 7 = 0. 
JUxV I 9s J 

Lemma 2.1. Under the above assumptions and (j2.6p - (|2.8j) . I(u,v) in (|2.ip /ias a <iua£ 
maximizing pair (u,v) G /C, where KL is the constraint set given in (|2.2p . 

Proof. The proof is inspired by [TJ [5]. Given any pair (u, v) G /C, we claim that I(u, v) does 
not decrease if v is replaced by 

(2.9) v*{y) = sup{s : <p(x , y , u(x) , s) < 0, \/x G [/}. 

In fact, by the continuity of (ft and m, for each y £ V there is some x G £/ such that 

4>(x,y,u(x),v*(y)) = > cf>(x,y,u(x),v(y)), 

since (u, u) G /C. By (|2.7p . u* > i>. Furthermore, (p(x , y , u(x) , v* (y)) < for all (x,y) G 
f/x V, so («,«*) G/C. 
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Since v* >v, by (|2.6p we have 

I(u,v*) >I(u,v). 
Similarly, if we define 

(2.10) u*(x) = sup{t : <f>(x,y,t,v*(y)) < 0, Vy G V}, 

then (u*,v*) G K, and 

I(u*,i;*) >/(«,«*) >I{u,v). 
Thus we do not decrease I(u,v) by replacing (it, v) by (u*,u*). The claim is proved. 

Define /Cc = K, n {u > Co}, where Co is a constant. The constant Co may be chosen 
negative and sufficiently small in the following context. We show that u* and v* are uni- 
formly bounded if (u,v) G K-c - Since u* > v,u > Co, by (J2-7)) we have for each y G V, 
s :=«*(y), 

4>(x, y, Co, s) < 0(x, y, u(a;), s) < 0, for all x G U. 

Then by (|2.7p again, there exists a constant Ci such that s < C\. This implies that 

(2.11) v < v* < d, 

we may choose Cl such that supy v* = C\. By a similar argument, there is another constant 
Co depending on <j> and C\ such that mfjju* = Cq. The constant Co > Co, since u* > u in 
U, and so («*,«*) G JC Co - 

We next deduce the lower bound of v* and the upper bound of u* by showing that u* 
and v* are locally Lipschitz functions. Consider two points in U, x\ ^ X2 and \x\ — x%\ < e 
is sufficiently small. There are two points yi, y2 G V such that 

cf)(x 1 ,y 1 ,u*(x 1 ),v*(y 1 )) = 0, 
<t>{x2,y2,u*{x2),v*{y2)) = 0. 

Then we have 

=4>{x 2 , V2,u* (x 2 ),v*(y 2 )) -4>{xi,y2,u*(x 1 ),v*(y 2 )) 

+ H x i^y2,u*(xi),v*(y 2 )) -cp{x 1 ,y 1 ,u*(x 1 ),v*(y 1 )) 

=4>t(x,y 2 ,u*,v*)(u*(x 2 ) -u*(xi)) - cj) x (x,y 2 ,u*,v*) ■ (x 2 - xi) 

+ 4>( x i^y2,u*{x 1 ),v*(y 2 )), 

where u* = 9u*{x\) + (1 — 6)u*(x 2 ),x = 9x\ + (1 — 9)x 2 , for some 8,6 G (0, 1). Noting that 

(f)(xi,y 2 ,u*(xi),v*(y 2 )) < 0, we have 

C 2 
u*{x 2 ) - u (xi) > -7r|x 2 - xi\, 

where the constants C 2 = sup(\d x (f>\ + l^^l), and C3 = minjinf d t (f>, inf d s (f)}. Due to p. 71) . 
the constant C3 > 5q is positive. On the other hand, replacing <p(xi,y 2 ,u*(xi),v*(y 2 )) by 
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4>(x2,yi,u*(x2),v*(yi)) in the above calculation, we have 

C 2 
u*(x 2 ) - u*(xi) < —\x 2 - Xl\. 

Therefore, the Lipschitz constant of u* on U is controlled by 

(2.12) \\u*\\Li P (u) < C 4 , 

where the constant C4 = C2/C3. A similar argument holds for v* as well, which implies 
that ||^*||i,ip(v) < C4. Hence, we have u* < Cq + C/±d\am.(U) and v* > C\ — CidiamiV) 
because of (|2.1ip . 

We conclude, therefore, that any pair (u, v) £ fCc may be replaced by a bounded, 

Lipschitz pair (u*,v*) G fCc without decreasing I. We now choose a sequence {(uk,Vk)} C 

/C(7 such that 

I(u k ,v k )^ sup I(u,v). 
(u,v)eKc 

By the above considerations we may assume that each (uk,Vk) is a bounded, uniformly 
Lipschitz pair, uniformly with respect to k, so there is a subsequence converging uniformly 
to a bounded, Lipschitz, maximizing pair (u,v) E ICco- 
Last, we show that when Co < is sufficiently small, 

sup I(u,v) = sup I(u,v), 

(u,v)£K.c (u,v)gK 

or equivalently, sup^ c / is independent of Co- By definition, one has sup^ c _ 1 1 > sup^ c /. 
So, it suffices to show the reverse inequality. Let (u,v) £ ICc -i be a maximizer such that 
I(u,v) = sup^ _j /, and {xk}k=i,--- ,N be a set of points in U. For a small constant e > 0, 

define 

u in U — U7v-B e (xfc), 

w + 2 in Uat B £ {xk)- 

Note that we may replace u by its mollification Uh = Ph * u, where ph is the standard 
mollifier function [3]. For simplicity, we assume u continuous in the sense that for h > 
sufficiently small, 

I(u h ,v) = I(u,v) + 0(Ne n ). 
Define 

v*(y) = sup{s : <f>(x, y, u(x), s) < 0, Vx G [/}, 
u*(x) = sup{t : <f>{x,y,t,v*{y))<0, Vt/eV}. 

Since the constraint function <fi is smooth and by (J2.7J) , except a set E C U and a set E' C V 
of measure |£| = |£'| = 0(iVe n ), 

u* = v- 2^ + 0(5) mV\E', 
u* = u + 2 + 0(6) mU\E, 
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where 5 := min^^,-{dist(xj,Xj)}. Therefore, by (|2.8p and the mean value theorem we have 
I(u*,v*) = I(u,v) + 2 [ f Ft - F &\ dl + (6)+0(Ne n ) 

J(U\E)x(V\E>) I 9s) 

>I(u,v)-C5-CNe n . 

As (u, v) E K,c -i, we may assume that inf[/ u = Cq — 1. Otherwise, one has mijj u = Co — tq 
for some constant tq < 1. This implies that sup^ c _ 1 1 = sup^ c _ t I, namely sup^- c I is 
independent of Co, and the proof is finished. By the definition, 5 will become small if the 
number of points N is sufficiently large so that we have (u*,v*) G K,c and 

supl>l(u*,v*)> sup I-C5-CNe n . 

Then, choosing e > sufficiently small we have 

sup I < sup I, 

by letting 5 — > 0, e — > 0, which implies that sup^ I is independent of Co, and the proof is 
finished. □ 

Remark 2.1. From the proof of Lemma 12.11 we conclude that there exist infinitely many 
maximizing pairs. In fact, if (u, v) is a maximizer and Co = infj/ u, then there is another 
maximizer in /Cc* +i, which is different from (u, v). 

Lemma 2.2. Let (u,v) £ /C be a dual maximizing pair in Lemma \2.1\ The equation 

(2.13) <f>(x, T(x),u(x),v(T(x))) = 

can be solved by a mapping T : U — )• V implicitly determined by the formula 

(2.14) 4> x (x, T(x),u, v) + 4>t(x, T(x),u, v)Du(x) = 0, 

at any differentiable point of u. Furthermore, if for any (x,y,t,s) E U x V x M?, 

(2.15) det [4> xy + 4> yt <&Du + o[ xs ® Dv + 4> ts Du ® Dv] / 0, 
the mapping T is uniquely determined by (|2.14p . 

The mapping T in Lemma 12.21 is called the optimal mapping associated to the dual 
maximizing pair (u, v). The inequality (J2.15P is a generalization of (A2) condition in optimal 
transportation [S]. 



Proof. Since u satisfies (|2.5p and v, 4> are continuous, for each x G £7, there exists some 
y =: T(x) G V such that 

(2.16) (f>(x,y,u(x),v(y)) = 0, 

(f)(x',y,u(x'),v(y)) < 0, 
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for any other x' G U. Let x G U be a differentiable point of u, by differentiation we have 

4> x (x, y, u, v) + 4> t (x, y, u, v)Du(x) = 0. 
If there exists another y ^ y in V such that 

4> x (x, y, u, v) + 4> t (x, y, u, v)Du{x) = 0. 

By the mean value theorem, 

(4>xy + 4>yt <8> Du + <j) xs (8) Dv + </>ts-D-u (8> -Du) ■ (y - y) = 0, 

where the matrix is valued at (x, y, u(x),v(y)) with y = ay + (1 — a)y for some a € (0, 1). 
This is a contradiction with (J2.15J) . since y ^ y. D 

Moreover, such an obtained optimal mapping T satisfies the following property, which is 
a kind of conservation of energy. 

Lemma 2.3. Let T be the optimal mapping associated to a dual maximizing pair (u,v). 
Assume the constraint function <p = 4>(x, y, t, s) is smooth and satisfies (|2.7p . Then for any 
h G C(V), there holds 

(2.17) 0=/ {-F t ^h{T(x)) + F s h{y)\d 1 . 

JUxV I Pt J 

Proof. Let h £ C(V) and |e| < 1 sufficiently small. Define 

(2.18) v e (y) = v(y) + eh(y) 
and 

(2.19) u e {x) = sup{t : <P{x,y,t,v e (y)) < 0, Vy G V}. 

Then (u e ,v e ) G K, and {uq,Vq) = (u,v). 

Since (u, v) satisfies (|2.5|) . by Lemma [2. 2 1 for every x G U the supremum (|2.5[) is attained 
at point j/o = T(x). We claim that at these points we have 

(2.20) uJx) - u(x) = -e^-h(T(x)) + o(e). 

To prove (l23(l . first we show that LHS < RHS. 
= (p(x,y ,u(x),v(y )) 
= 4>(x,yo,u{x),v e {y ) -eh(y )) 
= <P{x, 2/o, u{x),v e {y )) - e<p s h(y ) + o(e) 
= <j>(x, yo, u e (x) + u(x) - u e (x),v e (y )) - ecf) s h(y ) + o(e) 
= <^(a?,j/o,«e(a;),Ue(yo)) + <^t(u(x) - u e (x)) - e<p s h(y ) + o(e) 
< <pt(u(x) - u e {x)) - e<p s h(y ) + o(e). 
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By ([2771) we nave 

u e (x) - u{x) < -e-^h{y ) + o(e). 

To show LHS > RHS we use the fact that for any such x £ U there are points y e £ V 
such that the supremum in (|2.19p is attained. Thus 

> </>(x,y e ,u(x),v(y e )) 
= H x >ye,u(x),v e (ye) - eh(y e )) 
= <P{x,ye,u(x),v e (y e )) -ec/) s h(y e ) + o(e) 
= <f>(x, y e ,Ue(x) + u(x) - u e (x),v € (y e )) - e4> s h(y e ) + o(e) 
= <£(a:, y e ,u e (x),v e (y t )) + <fi t (u(x) - u t (x)) - e<f> s h(y e ) + o(e) 
= (f>t(u(x) - u e (x)) - e(f> s h(y e ) + o(e). 
Then by (|2.7p we have 

u e (x) - it(x) > -e-^-h{y e ) + o(e) 

<Pt 

= -e^% ) + e^ (%o) - h(y e )) + o(e). 

0t <Pt 

Since the supremum in (|2.5|> is attained at g/o> we have y e — )• yo as e — )• 0, and therefore, 
since h £ C(V), 

e- 1 (h(yo) - h{y e )) = o(e). 

This implies that LfT5 > ##S\ and (J272DJI follows. 

Next, since (it, i>) = (ito,^o) maximizes / among all pairs in /C, we obtain 

I(u e ,v e )-I(u,v) 

I) = lim 

e-^0 e 

/" F(x, y, it £ , i; £ ) - F(x, y, it, i;) 



c^O 



17 xV 






D 



2.2. Formulation of reflector problem. In order to formulate the near field reflector 
problem to an optimization problem, we need the notion of ellipsoid of revolution, which 
has a special reflection property: the light rays from one focus are always reflected to the 
other focus. 

In the polar coordinate system, an ellipsoid of revolution E = E(Y,p) with one focus 
at the origin, the other focus at Y, and focal parameter p £ (0, oo) can be represented as 



10 JIAKUN LIU 

E = {X Pe {X) : X G §"} by a radial function 

P 



(2-21) Pe(X) = — Y 

l-e(p)(X, m ) 

and 

(2-22) e(p) = \r^ _ M 

is the eccentricity, see [7] . Note that any such ellipsoid is uniquely determined by Y and p. 
If we regard p = piY) as a /oca/ function on $7*, we then have a family of ellipsoids. 

We recall that [6], for an admissible reflector T p , at each point Xp(X) G T p there exists 
a supporting ellipsoid, namely for some 7efl* 



(2-23) < ,„„ , '-«$<*.*->• 






In the following context, we also say p is admissible if T p is admissible. 
Next we define a set-valued mapping T p : — )• Q*. For any led, 

(2 24) Tp ^ = ^ Y E n * : Y is the f ° CUS ° f 

a supporting ellipsoid of T p at Xp(X)}. 

Note that at any differentiable point X of p, T p {X) is single valued and is exactly the 
reflection mapping. For any subset G C Q, we denote T p (G) = U XeG .T p (X). Therefore, we 
can define a measure \i# = /i Pi5 in such that for any Borel set G C ft, 



(2.25) p # (G) = / 5 di/. 

■/T P (G) 

Definition 2.2. An admissible function p is called a weak solution of the reflector problem 
if P>p,g = fdfi as measures, namely for any Borel set GcU, 



(2.26) fdfjL= gdv. 

JG JT P {G) 

The above definition was introduced in [6]. Obviously an admissible smooth solution is a 
weak solution, in that case, the reflector T p is naturally an envelope of a family of confocal 
ellipsoids of revolution. Therefore, the radial function p satisfies 

(2.27) p(X) = inf ^0. ^, led, 

Yen*l- e (p(Y)){X,jf ] y 

and for each Y G Ci* the ellipsoid Eyp(y) ls supporting to T p , we also have the focal function 
p satisfies 

Y 



(2.28) p(Y) = sup p(X) 



l-e(p(Y))(X, 



Ye n*. 
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Note that in (|2.27p for each led the infimum is achieved at some Y € £1* and in (|2.28p 
for each Y £ £1* the supremum is achieved at some X G VL. 

The relations ()2.27[ )— (j2.28|) between the radial and focal functions of a reflector T p are 
analogous to the classical relations between the radial and support functions for convex 
bodies, for example, see [TO]. Inspired by that and [II], we set r\ = 1/p. Then the pair (p, 77) 
satisfies the following dual relation 

(2.29) p{X) = inf , r , 

Y ^*r,(Y){l-e( V (Y))(X,^)) 

77(F) = inf 

X&Q 



where 77 is a Legendre type transform of p, |3J. 

Similarly to [14] , we can now formulate the reflector problem to a nonlinear optimization 
ipn ]) -ipr2 J) as follows. Set the functional 

I(P,V)= I F(X,Y,p, V ) 
Jnxn* 
(2.30) 



! f(X)\ogp + g{Y)(\ogr, + \og{l- — 

Jqxq* y i] i 



(X,Y) 



+ v / |y| 2 + 7 r 2 ' 

and the constraint set 

IC = {(p,r ] )eC + (n)xC + (n*) : <f>(X,Y,p,r})<0}, 

with the constraint function 

Y 



(2.31) <P(X,Y,p,ri) = log p + log V + logll-e( V (Y))(X^ 

In fact, by (|2.22p and 77 = 1/p it is easy to see that 

Y x (X,Y) 



(2.32) e( V (Y))(X, 



y\ v-' + vWW+v^ 



Lemma 2.4. Let (p,rj) be a dual maximizing pair of (|2.30|) f)2.31[) . and T be the associated 
optimal mapping. Then T = T p at any differentiable point of p, where T p is the reflection 
mapping in (|2.24p . 

Proof. We first introduce some geometric notation. By restricting to a subset we may 
assume that f2 is in the north hemisphere. Let X = (x, x n+ i) be a smooth parameterization 
of Oc S n , where x n+ \ = y/l — \x\ 2 and x = (x\, ■ ■ ■ ,x n ). 
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Denote di = d/dxi, e, = diX, and the metric gij = (e,, ej), where ( , ) is the inner product 
of R n+1 . By direct computations, for i,j,k,l = l,--- , n, 

(2.33) ej = 0, • • • , 1, • • • ,0, — l , 1 is in the ith coordinate, 

(2.34) g i:j = Sij + - _Vi 2 > y = kj ~ XiXj, where (^') = (gij)' 1 , 
and the Christoffel symbols are 

(2.35) I* = \g kl (d m + a^ fl - d m ) = x k Uj + y^i) • 
Denote e n +i = —X, the unit inner normal of § n at X. The Gauss formula is 



J 3 

where the second fundamental form 



(jj&i — 1 ij&k ~r 'T'ijS-n+li 



U.. _ X.. 1 ^J 

'41 U 21 ^ 



77 — u 27 ~ - I 19 

(2.36) l-\x\ 2 

= 9ij- 

Namely, one has that 

(2.37) djd = T%e k - 9lj X. 

The above equalities (I2.33p - (|2.37p can all be obtained by basic computations. 

Let p be a function defined on Q. The tangential gradient of p is defined by 

v. 

(2.38) Vp = ^ 9 ij e l d jP . 

Note that Vp(X) £ Tx§ n , the tangent space, so (Vp,X) = 0. By direct calculation 

(Vp,ei) = (g jk e j d k p,e i ) 
(2.39) 

= g 3 gijdkP = S ik d k p = dip, 

for all 1 < i < n. Let Dp = (dip, • • • , d n p) be the standard gradient of p. From (|2.33p . 
(12341) and (12381) . we have 

(2.40) Vp = (Dp,0)-(Dp-x)X, 

(2.41) | Vp| 2 = (Vp, Vp) = |Dp| 2 - (Dp • x) 2 . 

Let T p = {Xp(X) : X G fi} be the graph of p over Q. We claim that the unit normal of 
T p at Xp{X) is 

Vp - pX 

(2.42) 7 = - F H 

Vp 2 + IVpI 2 

Indeed, for z = 1, • • • , n, the tangential of T p at Xp{X) is 

ri = ^(xp(x)) = pei + (d iP )x. 
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From (|2.39p . for any i = 1, • • • ,n, the following holds: 

^ = / 2 !iV7 12 ^ + iPiP)X,V P - P X) 

Vp + l v H 

-.{pdip - pdip) = 0. 



Vp 2 + |vH 2 

It is obvious that |^y | = 1, thus 7 is the unit normal. 

At the differentiable point X of p, by (|2.24p . Y = T p (X) is the focus of the supporting 
ellipsoid of F p at Xp(X). Denote the reflected direction by Y r = \ Y ~_x P p \ • ^ f|2.42|) and the 
reflection law, 

Y r = X-2(X, 7 >7 
(2-43) 2pVp + (|Vp| 2 -p 2 )X 

|Vp| 2 + p 2 
Denote the length of reflected ray 

(2.44) d:=\Y-Xp\. 

Hence, we have 

Y = T p (X) = Xp + Y r d 

(2-45) 2pVp+(\Vp\ 2 -p 2 )X 



\Vp\ 2 + p 2 



d + Xp, 



and 

|Vp| 2 - p 2 
(2.46) {X , Y)=i \-A^L. + p . 

On the other hand, by differentiating the constraint function in (|2.31[) and the formula 
(|2.14p . we obtain 

P 1 - e(X, Y e ) 
where Y e = T(X)/\T(X)\, T(X) = Y is the optimal mapping. By noting that e* _L X, 
Vp _L X and (Vp, et) = d^p, from ()2.47p we have the decomposition 

(2.48) Y e = 1 ~ e{X ' Ye) Vp + (X, Y e )X. 

From (|2.13p . (|2.29[) and ()2.3ip . observe that at differentiable point X of p, there exists a 
unique supporting ellipsoid E of V p at Xp(X), with foci O, Y and eccentricity e. Note that 
the sum of length p = \Xp(X) — 0\ and length d = \Y — Xp{X)\ equals to the diameter of 
E, i.e. 

(2.49) p + d = diam(£). 

By the definition of eccentricity e, 

\Y\ \Y\ 



(2.50) 



diam(E') p + d 
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Combining (|2.48p and ()2.50p . one obtains the following equation for Y = T(X), 

(2.51) Y = P + d ~ {X ' Y) Vp + (X, Y)X. 

P 

It then suffices to show that Y = T P (X) in (f2T45|) is a solution of (J2~5Tj) . In fact, by (pT46|) 
we have 

T(X) = (d - d§^4) Vp/p + X P + Xd^ 



\V P \ 2 + p 2 J r/r r \Vp\ 2 + p 2 

: P Vp+(\Vp\ 
\Vp\ 2 4 
Xp + Y r d = T p (X). 



(2.52) 2pVp+(\Vp\ 2 -p 2 )X ^ 

~ p + WpTT? — d 



□ 



Proof of Theorem \1.1[ The proof essentially follows from [T3]. Let u = logp, v = logr]. 
In order to apply Lemma |2.1|, we need first to verify that F and 4> in (|2.30|) f)2.31[) satisfy 
the conditions (|2.6p - (|2.8p . For the constraint function <fi in (|2.3ip . it is easy to see that 
4> t = 1 > 0. By (^22]) and r? = 1/p, 




where 9 = (X, Y e ) G [—1, 1) due to (jl.5p . Since </> s is decreasing in 9, we have 

<P« > 1 , 

where the constant 9q < 1 depends on domains $7,0*. Set r = 1/(t/|Y|), e = y/l + r 2 
One has the second term in the above inequality 

< /. n o r =: ^W' 



^ 1 + r? 2|y|2 1_ e o 7r^72( 1 _^/f^72 +r ) 

where the function /i is decreasing in r. Thus, by the Taylor expansion of \J\ + r 2 near 
r = 0, 

/i(r) < lim h(r) = 1, for r > 0. 

r— >0 

Hence, we otain S > 5q and (|2.7[) holds, for a positive constant <5o- From (J2.3QH — (|2.32f) . one 

can see that Ft = 4>tf(X) and F s = <p s g(Y). Since / and g are both positive, we have the 
condition (12. 6p satisfied. The condition (12. 8p is an equivalent to the assumption (|l.ljl . 

Therefore, from Lemma 12.11 we obtain a dual maximizing pair (p, rj) £ /C in Theorem 
11.11 Then by the dual relation (|2.5p and (12.31H , one can see that p is admissible (J2.23P , and 
rj is the Legendre type transform of p as in (J2.29P . From Lemma 12.11 one knows that p is 
Lipschitz continuous. Actually, since an admissible function has supporting ellipsoid at any 
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point of its graph, it is semi-convex and twice differentiable almost everywhere [6] . Hence, 
by Lemma 12.41 T^ = T a.e., where T p is the mapping defined in (J2.24D , 

Next, we show that T satisfies the mearsure preserving condition (|2.26p . Since Ft = 
(j)tf(X) and F s = (j) s g(Y), by (12 .7p and applying Lemma 12.31 to T, we obtain that 



f(X)h(T(X))dfi = / g(Y)h(Y)du, 

n Jn* 

for arbitrary test functions h E C(Q*). Therefore, since T p = T a.e., we see that T p satisfies 
(I2.26p . namely p is a weak solution of the reflector problem. □ 



3. Derivation of equation 

We first derive the partial differential equation for the nonlinear optimization problem 
(|2.ip - (j2.2p in general. Let (u, v) be a dual maximizing pair of /. Assume that all the 
functions are smoothly differentiable at this stage. By a second differentiation of (|2.14j) we 
obtain 

=<j) xx + <j) X yDT + 2(j) x t ®Du+ ((f) xs ® Dv)DT 

+ (<t> yt ® Du)DT + 4> tt Du ®Du + (<f) ts Dv ® Du)DT + <f> t D 2 u, 

where each side is regarded as an n x n matrix valued at (x,y), y = T(x). 

Note that for every x E U, the equality (|2.13|) holds at point y = T(x), and for any other 
y' G V we have 

4>(x,y',u(x),v(y')) < 0, 
since (u,v) £ fC. Thus, at (x,T(x)) we have 



(3.1) 



By the assumption ([2] 
(3.2) 



Ty~^ 
> 0, we get 

Dv 



+ d>.Dv = 0. 



Combining (13. ip and (13. 2|) . we obtain the equation 

\(f) t D 2 u + 4> tt Du (g) Du + 2<fi xt ® Du + <fi xx \ 
= \4>xy + 4> xs ®Dv + 4> yt ®Du + cf> ts Dv (8) Du\ \DT\ 



(3.3) 



hence by ([2~7]) . 



f* 



y zy 



^ + 0j)t ® Du — —<t>y <8> Du 



\DT\ 



(3.4) 



det 



i 1 

D u -\ Du (g) Du H x t <g> -Dti + —4> xx 



det 



y zy 



+ 



V 



Du — t - </>3/ ® -Dw 



I det DTI 
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Equation ()3.4p is a second order fully nonlinear PDE of general Monge- Ampere type [3]. 
In the special case of optimal transportation (|2,3p - (|2.4p . equation (|3.4p becomes 



(3.5) 



Idet \D 2 u - D 2 c\ I = ldet.Di.cl P - . 

I L J I I - V ' p* O J 



For the derivation of the optimal transportation equation (|3.5p . see [9] for more. 

Using the notation from the proof of Lemma 12.41 we can now derive the PDE in the near 
field reflector problem by using the formula (|3.4p and constraint function (I2.3ip . 

Denote Y e = Y/\Y\, Dp = (dip,--- ,d n p) the gradient of p, and D 2 p = (didjp) the 
Hessian of p. By differentiating (|2.3ip . 



1 



1 



p p z 

e{ei,Y e ) 



>xt 



0, 



e(djei,Y e ) e 2 (e i) Y e )(e j ,Y e ) 



vx '~ 1 - e(X, Y e ) ' WXi *> " 1 - e(X, Y e ) (1 - e(X, Y e )f 

As in (|2.47p . at Y = T(X), where T is the optimal mapping in (|2.13p . we have 

(3.6) 

Therefore, 



dtp _ e(e i: Y e ) 
P 



l-e(X,Y e 
e(djei,Y e } 



(3-7) cj> XiXj - _ e{x ^ Ye ] 

and the LHS of equation (|3.4p becomes 



1 
? 



dipdjp, 



M(p) :-- 



(3.8) 



det 

det 



9 4>tt 2 1 

D 2 p + --Dp <g> Dp + —cj) xt ®Dp+—c 



d idjP - -d iP d jP - P l _ e{XtYe 



From ([OS]), 



(3.9) 



djei 



0,--- ,0, 



iL/^iL/j 



a/1 



(i-ix| 2 )^r 



.7.1 



In the special case Q,* C {y n +\ = 0}. 

Y 



ip ~ ' iyi — \*e,li > ie,m UJ. 



So, {djei,Y e ) = and (13.8P becomes 

M(p) = 



det 



D l p--Dp®Dp 
P 



Let u = 1//3. We have the standard Monge-Ampere operator as 

M(u) = |detD 2 u|. 
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In the general case when CI* is given by (II, 7|) . let us now calculate the term P i_Hxy) m 
(|3.8p . Let E be the supporting ellipsoid of T p at Xp(X), with foci 0,Y and eccentricity e. 
Recall that we have the relation ([2.501) . 

Therefore, 

e(djei,Y e ) pe{djei,Y) 



P 



(3.10) 



l-e(X,Y e ) \Y\-e(X,Y) 
p(djei,Y) 



from (13.9 



p + d-(X,Y) 
~P 



Vn+l 



$ii + 



1 — Ixl 



p + d- (X,Y) \x n+1 
Combining (I2.46P into (|3.10j) . we obtain 

,0-,-n pe(djei,Y e ) _ |V p\ 2 + p 2 ( y n +i\ Kr 

[ ] l-e(X,Y e )~ 2pd U+J '' 

where {Afij} is in (jl.lip . Actually, as one can see from (|2.36p . A/ij = gtj = hij is equal to 
the metric and the second fundamental form under the projection coordinates (|2.33p . 

Next, let us now calculate the length d = \Y — Xp\ appearing in (|3.1ip . Recall that 
Vp = g ij eidjp satisfies ([2^0]) - (p^T]) . Thus, from flgggp , we have 

d 



Vn+l 



(3.12) 

Therefore 
(3.13) 



\Vp\ 2 + p 2 

d 
|Vp| 2 + p 2 



(-2p(Dp ■ x)x n+ i + (|V/?| 2 - /) 2 )x n+ i) + px n+ i 
(\Dp\ 2 - (p + Dp- x) 2 ) x n+ i + px n+1 . 



d 



Vn+l 



P 



\Vp\ 2 + p 2 



x n+1 J \Dp\ 2 -(p + Dp-x) 2 ' 
Finally, combining (J3.13P into (|3.1ip we obtain 

pe(d jei ,Y e ) _ \Dp\ 2 -(p + Dp-x) 2 



(3.14) 



Vn+l 



Mi 



l-e(X,Y e ) 2p \px n+ i-y n+ i 

Using the notation (fl31)- (fTTill . a = \Dp\ 2 -{p + Dp- x) 2 and M = {Mij). We get the LHS 



of equation ([3? 



(3.15) 



M{p) 



det 
det 



D 2 p- -Dp® Dp 
P 

9 2 

D 2 p- -Dp® Dp 

P 



ay n +i 



2p(px n+ i - y n+ i 
o(l-t) 



-AA 



2tp 



-AA 



To compute the RHS of (|3.4p . one can directly differentiate the constraint function in 
(|2.3ip . but the computations are rather complicated. Instead, we recall a result in [6] in 
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the following: the Jacobian determinant of the reflection mapping T p is equal to 



(3.16) 



|det DT p \ = 2V ^n+ilV-01 



2 n p 2n+1 x n+1 \V^\ 



t n b/3 



a n+l 

t n b/3 



det 



D z p--Dp®Dp 
P 



a(l-t) 
2tp 



N 



,n+l 



M(p), 



where b,/3 are defined in (ll.9p -( [L"l~0~j) and ip is the defining function of £1* in (II. 7p . Alter- 
natively, one can obtain (|3.16p by differetiating the mapping T p in (I2.45p . 

On the other hand, from ([32 



(3.17) 



\det DT\ =M(p)<ft 



det 



1 



'xy 



-4> xs ®(f)y+ (f>y t <8> DU ■ 



'I* 



Du 



By Lemma El and Theorem O \det DT\ = |detLT p |. Thus 



^ T? 



det 



Jxy 



, + < 



'yt 



Du — c 



Du 



,n+l 



Note that we projected Q C S n on the n dimensional space {x\ 
codp, where dp is the surface area element of Q, uj = yl- 



t n bj3 

i 
2 . By Lemma 



2 n p 2n + 1 x n+1 \V^\' 

,x n ) in (|2.33|) . dx = 
andd22SD, 



(3.19) 



|det£>T| = |detL>T p 






Therefore, combining (|3.18[) f)3.19|) into (|3.4p . we obtain the equation 



(3.20) 



M(p) 



,n+l 



t"6/3 



/ 



2n p 2n+l w 2 5 | V ^| ' 

This completes the proof of Theorem 11.21 However, note that since we calculate the 
absolute value for the determinant, the matrix in M(p), (J3.15P has a different sign to that 
in [6]. 

Remark 3.1. Another special case of the reflector problem is the far field case [13] . Suppose 
a ray X is reflected off by T p to a direction Y. Set the functional and constraint function 
in (|230D - ([23T|) to be 



(3.21) 
(3.22) 



i(p,v) 

<j>(X,Y,p,r,) 



\ogp(X)f(X)+ / logr,(Y)g(Y), 

n Jn* 

logp + logri + log(l-(X,Y)). 



Similarly to Theorem 11.11 one can show that if (p, rj) is a dual maximizing pair of I, then 
p is a solution of the far field reflector problem. This formulation was previously obtained 
by Wang in [H] . 

The equation in the far field case can be directly obtained by using the formula (|3.4p and 
differentiating the constraint function (|3.22p . Here we remark that the far field equation is 
a limit case of (J3.15P for the near field one, [6]. 
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To see this, using our notations (|1.8p - (|1.10p . from (|3.13p we have the length of reflected 
ray d = \Y — Xp\ is equal to 

V Xn+i J \Dp\ 2 -(p + Dp-x) 2 'a 

Let's regard the target fi* = {rZ : Z G S7*}, where r is sufficiently large, and O^ is a 
domain in the south hemisphere of S n . In this case, the defining function in (|1.7p will be 
ip(Z) = r 2 — \Z\ 2 . Let g r be the light distribution on Q* under the same reflector T. Then 
when r is sufficiently large, r n g r — > g, and 

t\V?p\ -t a 



(3.24) /3|V^| 

(3.25) 



(Y-Xp) -Vijj d pb' 
r \Y\ d pb 

t t t a 

Sending r — > oo, from (|1.12p we obtain the equation for the far field case 



(3.26) 



det 



D 2 p- -Dp®Dp+-?-N 
p 2p 
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